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In this article, we present a study of high-energy neutrino emission in gravitational 
collapse. A compact star is treated as a complete degenerate Fermi gas of neutrons, 
protons and electrons. In gravitational collapse, its density reaches the thresholds 
for muon and pion productions, leading to high-energy neutrinos production. By 
using adiabatic approximation that macroscopic collapsing processes are much slower 
than microscopic processes of particle interactions, we adopt equilibrium equations 
of microscopic processes to obtain the number of neutrino productions. Assuming 
10% of variation in gravitational binding energy converted to the energy of produced 
neutrinos, we obtain fluxes of 10 MeV electron-neutrinos and GeV electron and muon 
neutrinos. In addition, we compute the ratio (< 1) of total muon neutrino number 
to the total electron neutrino number at the source and at the Earth considering 
neutrino oscillations. We approximately obtain the number of GeV antineutrino 
events (> 1) in an ordinary detector such as Kamiokande and total energy of neutrino 
flux (> 10 53 erg), as a function of collapsing star mass. 

PACS numbers: 97.60.Lf, 95.30.Tg, 95.85.Ry, 13.15. +g 



I. INTRODUCTION 



When the most of nuclear fuel of a star has been consumed, there is nothing that resists 
against the inward gravitational pressure and collapse occurs. At the end of gravitational 
collapse will born a compact star which can be white dwarf, neutron star or black hole. 
In white dwarf degenerate Fermi gas of electrons supplies the equilibrium pressure and in 
the neutron stars they are supported by pressure of degenerate Fermi gas of neutrons. The 
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collapse and supernova phase is thought to be accompanied by the most energetic neutrino 



burst CIS db neutron star or black hole is formed 



ft 



see reviews 



4J1. This picture has 

n 

been confirmed since the neutrino burst with mean energy of about lOMeV known as 
SN1987A, was detected by both the Kamiokande-II [5f and the 1MB 6| water Cerenkov 

nn 

detectors. In SN1987A a neutron star, whose mass is about 1.4 M , formed [7J, [8[. However, 
stars with the mass M > 3.2 M gravitationally collapse to black holes Recently 10], 
numerical simulations for the merger of binary neutron stars are performed incorporating 
a finite-temperature, Shen's, equation of state and neutrino cooling effect. They showed 
that for such system a hypermassive neutron star results and black hole is not promptly 
formed after the onset of the merger as long as the total mass of the system is smaller 
than 3.2 M Q . The neutrino luminosity of the hypermassive neutron star was shown to be 
~ (3 — 10) x 10 53 erg/s. As o ppo sed to the thermal lOMeV neutrinos, astrophysical objects 



like supernova core 



collapses 
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13[, magnetars [14J, neutron stars 



supernova remnants 

and pulsars [15|, ll6| can produce high-energy (> 1 GeV) neutrinos, see reviews (I7l4l9l|. The 
production mechanism of high-energy neutrinos is the decay chain of charged pions and 



muons 



7T -)•// + r/ M e + r/ M + r/ M + v e , 
7T + n + + ->■ e + 



(1) 



these pions are produced by accelerated protons 20|, |21[ to high energies in the system via 



the photon-hadron interactions 



22| and inelastic hadronic processes [23] 



p + 7 -)■ A + -^n + 7r" 
p + p -> p + n + 7r + , 

where 7 is the target photon. The most dominant channel for pion production in Eq. 
is near the A resonance 



(2) 
(3) 



221 ] . Kilometer-scale neutrino detectors such as IceCube are con- 



structed to observe these high-energy neutrinos 

In this work we study another source of high-energy neutrinos from gravitationally col- 
lapsing stars to black holes. We compute number and mean energy of produced neutrinos in 
collapsing stars to black holes. To compute these quantities one has to solve rate equations 
which incorporate both, macroscopic hydrodynamic gravitational collapse process effects and 
microscopic particle interactions which produce neutrinos. However, time scale of particle 
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interactions is much shorter than time scale of collapsing processes. As a result, gravita- 
tional collapse processes can be considered adiabatic in comparison with particle interaction 
processes. Based on these argumentations, we divide collapse process into infinitesimal steps 
where in each step the system is in its thermodynamic equilibrium state, which we call the 
adiabatic approximation 25[. In addition to these thermodynamic equilibrium equations 
there are conserved quantities, which are total electric charge number, total baryon number 



and total lepton number of each flavor 



261 ] . We consider the collapsing star clS db complete 



degenerate Fermi gas of neutrons, pro tons and electrons because the temperature is much 
smaller than their Fermi energies 0, Solving equilibrium equations analytically, we 
obtain the number of each spices of particles as a function of collapsing star density. When 
collapse proceeds density of star reaches to the amount that muon then pion can be produced 
in the system. As a result, in addition to the electron neutrinos which are already produced 
by /3-processes, in companion to the muon production, muon neutrinos are produced. Using 
number of produced neutrinos and assumption of 10% 28|] of gravitational binding energy 
converted to the energy of produced neutrinos, we compute the mean energy of neutrinos. 
We study neutrino oscillation effects on the primary flux, and compute the number of an- 
tineutrino events in an ordinary detector such as Kamiokande. The total energy of neutrino 
flux is another relevant quantity to observations and we compute it in this article. 

The paper is organized as follows: In Sec. |TT1 we solve equilibrium equations and obtain 
number of each species of particles as a function of star density. In Sec. IHIl number of 
produced neutrinos is computed. In Sec. IIVI we compute mean energy of produced neutrinos. 
In Sec. |V] we compute the ratio of total muon neutrinos to the total electron neutrinos at 
the source and at the Earth, considering neutrino oscillation effects; and the number of 
antineutrino events in an ordinary detector such as Kamiokande obtained, as well as total 
energy of neutrino flux. Finally in the last section we summery and discuss the results. 



II. EQUILIBRIUM STATES 

In this paper, we approximately consider compact stars, white dwarf and neutron stars, 
as a complete degenerate Fermi gas of neutrons, protons and electrons for the reason that the 
temperature is much smaller than their Fermi energies. In equilibrium state, each reaction 
between particles is balanced by its inverse reaction and leads to the chemical equilibrium 



condition [26|, U_ J 

VidYi = 0, (4) 

i=particlcs 

where \i\ is the chemical potential of the particle-i, and Yi = n^/ns is the concentration of 
the i-species of particle with respect to the baryon number density ns such that rij is the 
number density of the particle-i. Therefore in equilibrium state the concentrations of the 
particles are not independent of each other, they are governed by the chemical equilibrium 
condition (J3J). In addition, total electric charge number, total baryon number A and total 
lepton number of each flavor, must be conserved. Suppose that high-energy neutrinos are 
allowed to escape from the system, neutrino number densities and consequently their Fermi 
energies can be approximately set to zero in Eq. fl4]) for equilibrium state. Equation (j4j), total 
charge and baryon number conservations, are sufficiently determine all numbers of particles 
in equilibrium states. 

A. /3-processes 

In the density region above the white dwarf density p w d = 8.73 x 10 -8 p nuc , where the 
nuclear density is p nnc = 2.8 x 10 14 grcm~ 3 , electrons are ultra relativistic and have enough 
energy to overcome the threshold 

E F , C > {m n - m p ), (5) 

then inverse /3-decay reaction 

e + p -> n + z/ e , (6) 

can proceed. In equilibrium state, the inverse /3-decay reaction is balanced by its back 
reaction i.e., /3-decay reaction 

n -> e + p + z/ e . (7) 

Variations of particle numbers due to the reactions and ([7|) respectively are given by 

dN Vc = dN n = -dN p = -dN e , (8) 
dN Pc = dN p = dN c = -dN n , (9) 
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where iV,-: 



n, p,e, v e and z/ e , respectively are total numbers of neutrons, protons, elec- 



trons, electron-neutrinos and electron-antineutrinos. Thus the chemical equilibrium condi- 
tion (jlj) requires following relation between chemical potentials of particles 

/in = Pp + fie- (10) 

For a complete degenerate system at zero temperature T, chemical potentials pi = 
mf + p Fi , and Fermi momentum p F ^ are known as functions of particle number densi- 
ties 71;, 



'ATA 1 / 3 



-1/2 



.4 



ft 



1/3 



(11) 



where we adopt the approximation for spatial homogeneity and the numerical coefficient 
£ = 0.173. Total charge neutrality and conservation of total baryon number respectively are 
given by 



N p = N e , 
A = N Q + N p . 

We obtain the solutions N p , N c and iV n to Eqs. (JT2D and (113"]) 



1 + ^(m 7r ) ( 



-2/3 



-3/2 



i + e 



2 2 2 

— nip — m c 
ml 



-2/3 N 



-4^ 



2 / m p m c 



mi 



N e (p) = N p (p), 
N D (p) = A-N p (p), 



(12) 
(13) 



-4/3 



3/2 



(14) 
(15) 
(16) 



such that in obtaining Eq. (1141) . only leading order terms i.e., (N p /A) 2 ^ 3 are taking into 
account and higher order terms are neglected. Similar formulas to Solutions (|T4l)-(fT6l) can 



be found in 



26 



27| . Solutions (fT^l) -(|T6l) completely determine the numbers of protons, 



electrons and neutrons as a function of the density p of compact stars. In Fig. [IJ we plot the 
ratio A/N p as a function of the density p, which is an important quantity to show charged 
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FIG. 1: The quantity A/N p is plotted as a function of the stellar density p in unit of nuclear 
density p nuc - 

compositions of compact stars at the density p. As shown in this figure, there are two 
main phases for neutrino productions by /3-processes. In the first phase, A/N p increases, 
through the inverse /3-decay reaction ([6]) electrons and protons having been converted into 
neutrons and electron- neutrinos, the later escape from compact star. The ratio A/N p reaches 
a maximum when the density reaches to the value 



p r ~ 21.75 /0b 



4((m n — m p ) 2 — ml) 



mi 



3/4 



2.8 x 10~> r 



(17) 



which is called the transition density [26|, |27j. In the second phase beyond the transition 
density (I17j) the ratio A/N p decreases monotonically, through /3-decay reaction (JTj) neutrons 
having been converted into protons, electrons and electron-antineutrinos escaping from com- 
pact star. We assume this second phase to end at the density at which other processes 
become important. 



B. Production of muons 



When the stellar density p reaches to the value p^ = 2.7 p nuc , the Fermi energy of electrons 
exceeds the muon mass, i.e., Ep c > m M , then it is energetically favorable for electrons to 
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turns into muons via 



And muons can decay to electrons 



e + u e -> (i + v„. 



H -> e + i/„ + v c 



(18) 



(19) 



so that muons and electrons are in equilibrium 



26] . The variations of particle numbers due 



to the reactions (jlsjl and ( 1191) respectively are given by 



dN- Vi 



dN„ 



-dN D 



-dN P 



dN„„ = dN 



-dN, 



(20) 
(21) 



where N; i — fi, and z/ M , respectively are total number of muons, muon-neutrinos 
and muon-antineutrinos. The chemical equilibrium condition (jlj) gives the /3-equilibrium 
condition (flOj) and the following relation between chemical potentials of electrons and muons 



26] 



/i e = AV (22) 
The conservations of total baryon numbers ffl3|) and total charge neutrality lead to 

N P = N C + TV (23) 



In terms of the particle densities, Eqs. (jTOj) . (jTBj) . (1221) and ( 1231) form a sufficient set of 
equations to solve for finding particle numbers as a function of the stellar density. We 
obtain the solutions for the numbers of electrons, muons, protons and neutrons as functions 
of the stellar density p 
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N n (p) = A - iVp(p), 



(24) 

(25) 

(26) 
(27) 
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FIG. 2: The quantity A/N p is plotted as a function of the density in the region p> = 2.1 p nVLC 
in solid line. Dashed line indicates the result obtained from /3-processes. 

such that in obtaining Eq. (I!H|) . only leading order terms i.e., (N c /A) 2 ^ 3 are taking into 
account and the muon number are neglected in compared to the electron number N e . 
We plot the ratio A/N p in Fig. [2], and find it decreases faster than what is predicted by 
/3-processes due to muon production. 




C. Production of pions 

If one ignores the effects of the strong interactions between pions and nucleons, the 
criterion for the formation of negatively charged pions in a compact star via the reaction 

n^p + TT, (28) 



is given by 26] 



A*n - A* P = /V (29) 

At low temperature T ~ 0, produced pions have approximately zero kinetic energy thus 
their chemical potential is close to the pion rest mass fi n ~ m n , that determines the pion 
production threshold at the density p n = 4.6p nuc . In equilibrium sate, produced pion can 
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decay to muon which balanced by its back reaction 



7T 



H + v 



i-i- 



(30) 



hence chemical equilibrium condition (J4j) requires 



(31) 



The chemical equilibrium condition leads to Eqs. ( )29|) and (I3ip . which are accompanied by 
the /3-equilibrium f lTUj) and electron and muon equilibrium (122]) . As a usual in addition to 
the baryon number conservation ( JTBl . charge neutrality leads to 



iV p = iV c + N„ + N n , 



(32) 



where N n is the number of produced pions (tt~)- By using expressions of chemical potentials 
in terms of particle densities, and solving Eqs. ( ITUj) . ( ITBl . ( |22|) . ( 129|) . ( 1311) and ( 1321) . we obtain 
solutions for the numbers of electrons, muons, protons, neutrons and pions as a function of 
the stellar density p 



3/2 (_P_ 
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(33) 
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N p (p) ~ —A 
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iV n (p) = A-N v (p), 

NM = iV p (p)-iV e (p)-7V», 



(35) 

(36) 
(37) 



such that in obtaining Eq. fl35l) . only leading order terms i.e., (iVp/A) 2 / 3 are taking into 
account and higher order terms are neglected. In Fig. [HI we plot the radio A/N p in this 
region p > p n , and find the production of negatively charged pions causes a faster decreasing 
of this ratio, comparing to the result obtained from /3-processes and muon production. 
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FIG. 3: The quantity A/N p is plotted as a function of density in the region p > p w = 4.6p nuc in 
solid line. Dotted and dashed lines respectively show what are predicted by /^-processes and muon 
production. 

III. NEUTRINO PRODUCTION IN GRAVITATIONAL COLLAPSES 



Stellar core density increases during gravitational collapses. As studied in the Sec. HTIpar- 
ticle numbers JV niP]Cift7r are functions of the stellar core density p, thus change in collapses. 
These variations of particle numbers are due to the /3-processes (EJ), (E3) and processes of 
muon-production ( fT8l) . ( TT9l) and pion-production ( |28i) . ( )30i) . Because the time-scales of these 
microscopic processes are much shorter than the time-scale of macroscopic processes, such as 
gravitational collapse and hydrodynamical processes, we approximately treat gravitational 
collapse and other macroscopic processes as adiabatic processes. Namely, gravitational col- 
lapse process can be divided into infinitesimal steps and at each step the system can be 
considered in its equilibrium state. Hence the results of the Sec. [Ill are applicable at each 
step of gravitational collapse. 

Let M be the mass of collapsing stars. When M > 3.2 M@, black holes are formed at the 



Schwarzschild radius FL = 2GM of the star 



where G is gravitational constant. This 



indicates a maximal density p s = M/(^wR^), which is plotted as a function of the mass M 
of collapsing stellar cores in Fig. HJ to compare with the density thresholds p^ and p w of 
muon and pion productions discussed in the Sec. [TlJ As a result, we conclude the following 
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FIG. 4: The density p at the Schwarzschild radius of the star p s , is plotted as a function of collapsing 
star mass M, in unit of Solar mass Mq, in solid line. Dashed and dotted lines respectively indicate 
the muon and pion production threshold densities p^ and p n . 

three cases 

1. 3.2 M & < M < 3.8 M Q , the density thresholds p M and p n of both muon and pion 
productions can be reached; 

2. 3.9 M Q < M < 4.9 M , only density threshold p M of muon production can be reached; 

3. M > 5 M & the density thresholds p M and p n of both muon and pion productions cannot 
be reached. 

According to these three cases and the density thresholds discussed in the Sec. [TTl we are 
ready to calculate the numbers of neutrino produced in different phases. 

A. Phase a 

In the phase a from white dwarf density p w d to the transition density pt HTTT) . Eq. (JH|) 
gives dN Uc = —dN c , thus the number of produced electron-neutrinos is 

K c = - dN c = A c (p wd ) - A c (p T ) ~ A/2, (38) 
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FIG. 5: The number of electron-antineutrinos Np a , in unit of Solar baryon number Aq ~ 1.19x 10 57 , 
produced in the phase b, is plotted as a function of collapsing star mass M. 



where numbers of electrons, protons and neutrons are 



iVc(Pwd) = iV p (p wd ) ~ iV n (p wd ) ~ A/2, 
at the white dwarf density p wd , and Fig. [1] shows 



(39) 



N c (p T ) = N p (p T ) ~ 1.34 x 10~ 4 A, 



(40) 



at the transition density px- 



B. Phase b 



The phase b of neutrino production begins from the transition density px to the density 
threshold p M for cases 1 and 2, and to the maximal density p s for case 3. The number of 
produced electron-antineutrinos, due to the /3-decay reaction (JJJ), is equal to the variation of 
electron number dNp c = dN e shown in Eq. Q). In cases 1 and 2, we obtain the total number 
of produced electron-antineutrinos 



No 



/ dN c = iV e (p M ) - 7V c (px) ~ 1.13 x IQ~ 2 A, 

J Pt 



(41) 
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FIG. 6: The number of electron-antineutrinos Np c and muon-antineutrinos produced in the 
phase c are plotted as a function of collapsing star mass M, respectively in solid and dashed lines. 



where Nd^p^) is given by Eq. ( 1241) . In case 3, we obtain 



dN e = N e (p s ) - N e (p T ), 



(42) 



Pt 



where N e (p s ) is given by Eq. ( fl~5l) . In Fig. El we plot Eqs. (I4"T1) and ( H2l) . which shows 
increases for cases 1 and 2, and decreases as collapsing mass M increases for case 3. 

C. Phase c 



The phase c of neutrino production begins from muon production density p^ to the density 
threshold p^ for cases 1, and to the maximal density p s for case 2. In this phase in addition 
to the /3-processes, reactions ffl~8]) and f|T9|) contribute to neutrino productions. Total number 
variations of electrons and muons, given by Eqs. ( 124)) and ( )25|) . are equal to the number of 
produced electron and muon antineutrinos, i.e., dN Pc = dN e and dN Pfi = dN^, due to the 
conservations of electron and muon flavor numbers. In case 1, we obtain the numbers of 
produced electron-antineutrinos N Pc and muon-antineutrinos iVp 



Nr, 



P» 
pit 



dN e = N C ( P7T ) - N c ( Pfl ) ~ 5.7 x 10" 3 A, 



dN^ = N^(p n ) - ~ 5.4 x 10~ 3 A 



(43) 
(44) 
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FIG. 7: The number of electron-neutrinos iV„ e and muon-neutrinos N v produced in the phase d 
are plotted as a function of collapsing star mass M, respectively in solid and dashed lines. 

where we use Eqs. ( }24"1) and (1251) . In case 2, we obtain the number of produced electron and 
muon antineutrinos 

N Pc = dN c = N c (p s ) - N c ( Pfl ), (45) 

Nv,= P dN„ = NM-NM, (46) 

where N e (p s ) and A^(p s ) are given by Eqs. ( 124]) and ( |25|) . The results ()43 i) -( l46|) are plotted in 
Fig. [6l which shows the numbers N Dq and A^p M increase for case 1 and decrease as collapsing 
star mass M increases for case 2. 

D. Phase d 

The Phase d of neutrino production begins from the threshold density p n for pion produc- 
tion and ends to the maximal density p s . This phase can only be reached in case 1. Equations 
( |33l) and ( 1341) show that by increasing the density p, electron number N e and muon num- 
ber Nfj, decrease. Therefore electron and muon neutrinos are produced, dN Uc = —dN c and 
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FIG. 8: Mean energy E Uc , in unit of MeV, of produced electron-neutrinos in the phase a is plotted 
as a function of collapsing star mass M. 



dN v = —dNp, and the total numbers are 



N„ 



Ms 



dN c = N c (p n ) - N e (p s 



Pit 

Pb 



dN^ = NuM-Npips), 



(47) 
(48) 



Pit 



which are plotted by Fig. [7J we find the number of neutrino productions are decreasing 
functions of collapsing star mass M. 

IV. MEAN ENERGY OF NEUTRINOS 



The gravitational energy of stellar core with mass M and radius R is given by 



R 

I 4vrr 2 
'o 



1 



2GM(r] 



-1/2 



p(r)dr, 



(49) 



where M.(r) = 4irr' 2 p(r')dr' [27]] . Assuming uniform density p = M/^irR 3 ), we have 

3 



U = M 



1 — —a 3 ^arcsin(a) — aVl — 



2GM 
R 



During collapses from radius R to R — dR, the variation of gravitational energy 



dU 



3M 1 
yT^ 



1 + 3a 3 ( yl — a 2 arcsin(a) — a 



dR, 



(50) 



(51) 
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FIG. 9: Mean energy Ep c , in unit of GeV, of produced electron-antineutrinos in the phase b, is 
plotted function of collapsing star mass M. 



is converted to the kinetic energy and internal energy of the particles. 

To compute the mean energy E v of emitted neutrinos, we assume 10% o 
of gravitational energy dU converted to the energy of produced neutrinos 



total variation 
. Suppose that 



N u neutrinos are emitted in gravitational collapse process from the radius R\ to the radius 
i?f, observed mean energy of neutrinos can be calculated by 

0.1 ft* 



where Z is the gravitational red shift factor 



Z(R)dU, 



(52) 



Z(R) 



1 



2GM 
R 



(53) 



Using Eqs. (!5T!) - (l53l and the number of produced neutrinos in each phases discussed in the 
Sec. IIIIl we obtain mean energies of neutrinos. In the phase a, the mean energy of produced 
N Uc , see Eq. (1381) . electron-neutrinos is given by 

0.1 rRT 



E, 



Z(R)dU, 



(54) 



where i? w d is the radius at which the collapsing core density p is equal to the white dwarf 
density p w d, and Rt is the radius at which the density reaches to the transition density px, 
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FIG. 10: Mean energy electron-antineutrinos and muon-antineutrinos Ep ct p fi , which are produced 
in the phase c, is plotted as a function of collapsing star mass M. 

see Eq. (1 17p . We plot this result (l54"j) in Fig. [Sj and find the neutrino mean energy is about 
10 MeV, which is in agreement with neutrino energy observed in supernova explosions. In 
the phase b, the mean energy of produced electron-antineutrinos for cases 1 and 2 is given 
by 

0.1 [*» 



E, 



Z(R)dU, 



(55) 



where N Pa is calculated by Eq. ( 141 j) and is a radius at which the density p reaches to p^. 
Whereas for case 3, Eq. ( 153]) becomes 

0.1 rB * 



Z(R)dU, 



(56) 



"vc J Rt 

where N Pc is calculated by Eq. (1421) . The results (|55|) and (|56|) are shown in Fig. [9j Compared 
with results in the first phase, the electron-antineutrino mean energy about GeV is larger 
than the electron-neutrino mean energy, because in this phase the variation of gravitational 
energy is larger and the number of produced neutrino is smaller, see Fig. [5] and Eq. (138|) . 
In the phase c, the mean energy of produced electron-antineutrinos and muon-antineutrinos 
for case 1 is given by 

0.1 rR * 



E- - 



Z(R)dU, 



(57) 
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FIG. 11: Mean energy electron-neutrinos and muon-neutrinos E UcjI/)1 , which are produced in the 
phase d, is plotted as a function of collapsing star mass M. 

where R n is the radius at which the density p reaches p n , neutrino numbers N Pc and N D are 
respectively given by Eqs. (jT3j) and (jHJ). Whereas for case 2 Eq. floTl) becomes 



E 



0.1 



M *"" (N Pe + N p , jRu 



Z{R)dU, 



(58) 



where neutrino numbers N& e and Np are respectively given by Eqs. (1451) and (|46|) . The 
results (1571) and (1581) are shown in Fig. [UJ In the phase d, the mean energy of produced 
electron-neutrinos and muon-neutrinos is given by 

0.1 



E„ 



Z(R)dU, 



(59) 



(N„ e + N UfijjRiv 

where the neutrino numbers N Ue and N v are respectively given by Eqs. (ITTj) and (j48|) . The 
result (159|) is shown in Fig. [TTJ We find that the mean energy neutrinos produced via muon 
and pion productions is about GeV, see Figs. [10] and [TTJ We also calculate the mean energy 
of total produced neutrinos throughout phases a, b, c and d 

0.1 rRs 



N„ 



Z(R)dU; 



Vc ISc. I^c i^jj, Vq. l/u ' 



(60) 



where the superscripts show the related neutrino production phases. The result (1601) is 
shown in Fig. [121 an d we find neutrino mean energy is about 90 MeV, which is larger than 



19 



93. 



92.5 

> 




91.5 



91. / 

5 10 15 20 

M/M Q 

FIG. 12: Mean energy E Ut p of total produced neutrinos and antineutrinos in the phases a, b, c and 
d, is plotted function of collapsing star mass M. 

neutrino mean energy lOMeV in supernova explosion. This is due to more gravitational 
energy converted to neutrino energy when stars gravitationally collapse to black holes. 



V. OSCILLATION AND DETECTION OF NEUTRINOS 



There are several observations, see [30j and references there in, indicate that neutrino has 
mass. Thus a neutrino with flavor £ is a superposition of the physical fields v a with different 

masses m n 



We 



> = X^K>> ( 61 ) 



where U is unitary matrix which signifies neutrino mixing [31]. By considering a simple- 
minded approach to the propagation of this state, one can assume that the three-momentum 
p of the different components in the neutrino beam are the same. However, since their 
masses are different, the energies of all these components cannot be equal, and given by 



E a — ^/p 2 + m a- I n t ne situation that we are studying, neutrinos are extremely relativistic, 
so that we can approximate the energy-momentum relation as E a ~ |p| + m 2 a /(2\p\). After 
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time t, or distance L, the evolution of the initial beam of Eq. (IBTj) gives 



tofi 



iE Q t 



u a ). 



(62) 



Thus the probability of finding a vg in the original V(_ beam, at any distance L, is 



P, 



\(MML))\ 2 

6 U > - 4 ^ u eaUe> a U e(3 U il( 3 sin 2 

a>/3 



Am 



q 3 



AE 



(63) 



where we use the notation \p\ — E and 



Am 2 a/3 = 



2 2 



(64) 



n writing Eq. (J63l) CP-conservation has been assumed when the mixing matrix [/ is real 



3l| . For three generations, the most general mixing matrix ignoring CP violation is given 
by 



U 



( 



C12C13 — S12C13 S13 

S12C23 + Ci2'S23'Sl3 C12C23 — S12S23S13 — S23C13 
S12S23 - C12C23S13 C12S23 + S12C23S13 C23C13 



7 



(65) 



where c^- = cos% and Sij = sin^. Analysis of the most recent data 32J on the oscillation 
parameters indicate that 



sin 2 (20 12 ) = 0.87 ±0.03, 
sin 2 (20 23 ) > 0.92, 
sin 2 (20 13 ) < 0.15, 



Am 2 ! = (7.59 ± 0.20) x 10" 5 eV 2 , 
Am 2 2 = (2.43 ± 0.13) x 10~ 3 eV 2 , 



(66) 



using these data, we chose 



0i 2 = 9 ~ 34.5°, 023 = 45°, 9 
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0. 



(67) 



For the mass squared differences Am 2 , which are given in Eq. fl66l) . energy dependent terms 
in the probability expression (1631) . even for high-energy neutrinos (E v > 1 GeV) traveling on 
cosmologically scales L ~ Mpc, are very fast oscillating functions and can be averaged to 



sin 



Am 2 \ ZlGeV 



IMpc J 



1 

^2' 



(68) 
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FIG. 13: The quantities <f> and cf)^' are plotted as a function of collapsing star mass M, respectively 
in solid and dashed lines. 



On the basis of results of the Sec. IHIl it is clear that primary flux of 10 MeV neutrinos 
contains only electron-neutrinos, and their numbers given by Eq. fl38|) ; and primary flux of 
GeV neutrinos contains electron and muon neutrinos, and their total numbers are given by 



(c) 



Ma 



N, 



(c) 



(69) 
(70) 



where these numbers can be read from Eqs. (|4ip - (j48|) . However oscillation change these num- 
bers, namely produce tau neutrinos which are absent in the primary flux. Using Eqs. (163]) - 
(!68|) after considering neutrino oscillations, the numbers (|69l and (1701 become 



AC (D) 



[1 - ^ sm 2 (26))Af c + l - sin 2 (2^)A^, 



i(l - ~ sin 2 (2fl))A/; + ~ sin 2 (2^)A4. 



(71) 



(72) 



Because water Cherenkov detectors, such as Kamiokande, can distinguish electrons from 
muons but cannot distinguish particles from anti-particles 30], it is useful to investigate, 
the ratio of the total produced muon neutrinos Af u to the total produced electron neutrinos 



22 



10.0 



5.0 - 



2.0 - 



1.0 



0.5 - 



0.2 - 
0.1 - 



10 



15 



20 



M/M Q 



FIG. 14: In a gravitational collapse process, the number of GeV antineutrino events which can be 
observed by Kamiokande detector is plotted as a function of collapsing star mass M. Thick and 
Thin lines respectively show numbers of electron and muon antineutrino events without considering 
oscillations. Dashed line shows numbers of electron-antineutrino events by considering oscillations. 
Dotted line shows number of muon and tau antineutrino events by considering oscillations. 



AC i.e., 



K 
A/1 



(73) 



Using Eqs. ( 1711 and (172]) after considering neutrino oscillations the quantity becomes 

, (D) _ A/f } _ |(l-|sin 2 (2g))0+lsin 2 (2g) 

9 " A; (D) l-±sin 2 (2£) + ±sin 2 (2£)0 ' {< ) 

which shows if <p — 2, regardless value of mixing angle 6, flux becomes homogenous i.e., 
</>( D ) = 1. For the flux of lOMeV neutrinos and collapsing stars with the mass M > 5M , 
we have (p = as a result = 0.39. Fig. [13] shows <p an d (f>^ for flux of GeV neutrinos, 
and Fig. [14] shows antineutrino events that can be detected on the Earth. These results can 
be relevant to observations. 
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A. Detection of neutrinos 

If one is interested only in the charge current processes &e+p — > £ + +n and v^+n — > £~+p, 
so that the neutrino is converted to a charged lepton that can be detected. From the 
experimental view point, only the reaction vi + p — > £ + + n is interest, because it is not 



possible to build a target containing enough free neutrons, that would anyway decay [30 ] . 
Enough free protons are obtained using targets made of water H 2 such as the Kamiokande 
detector. For this reason, analogously to Ref. [loj], we consider only antineutrino events 
without taking into account the numbers of lOMeV electron-neutrinos, which are produced 
in the phase a, and GeV electron-neutrinos and muon-neutrinos which are produced in the 
phase d. The number of observed events in a detector depends on the number of emitted 
neutrinos N u , cross section of neutrino-nucleon interaction a, number of detector target 



nucleons N and distance from collapsing stars to the Earth L. Thus it is determined as 31] 



N cvcnt = ^°N. (75) 

Here we compute number of GeV antineutrino events, because we are interested in high- 
energy neutrino detection. The cross section of GeV antineutrinos interaction with a proton 



approximately is a ~ 10 37 cm 2 30j]. The Kamiokande detector contains 2.2KTon of water 



i.e., N ~ 1.5 x 10 32 hydrogen nucleus or proton 3jJ. For a collapsing star at distance 
of L = 1 Mpc from the Earth, the number of observed GeV antineutrino events in the 
Kamiokande detector iV cvcnt > 1, as shown in Fig. [TU depending on collapsing star mass. 

Another important information is the total energy of neutrino flux. Total energy of 
10 MeV electron-neutrinos flux produced in phase a is 

e Me v = N Uc x E Ua , (76) 

where N Uc and E Ve are respectively given by Eq. ( )38|) and shown in Fig. |HJ Total energy of 
GeV electron and muon neutrinos fluxes produced in phases b, c and d is 

ecev = x Ef c + (Af + N$) x E%» + « d) + N$) x E^, (77) 

where , , , N^f , are shown in Figs. (MH and 4? , > in Fi S s - ™ 

Then, the total energy of neutrino flux e is the sum of eMev (EHJ) and ecov (ED)- Eqs. (1761) . 
( 1771) and total energy e of neutrino flux are presented in Fig. [151 which shows the total energy 
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FIG. 15: In unit of 10 53 erg, the total energy of neutrino flux e (solid line), the total energy of 
lOMeV electron- neutrino flux eMeV (dashed line), and the total energy of GeV electron and muon 
neutrino flux eQeV (dotted line) are plotted as a function of collapsing star mass M. 

of emitted neutrinos e > 10 53 erg, depending on collapsing star mass. The total energy of 
flux presented in Fig. [15] can be relevant for observations. 




VI. SUMMARY AND REMARKS 



In this article, we present a preliminary study of neutrino productions in gravitational 
collapse of stellar cores. For this propose we consider the collapsing star as a completely de- 
generate Fermi gas of neutrons, protons and electrons, because temperature is much smaller 
than their Fermi energies. To compute neutrino numbers and mean energies one has to solve 
rate equations, which incorporate both macroscopic hydrodynamic collapsing processes and 
microscopic particle interactions. However, time scale of macroscopic collapsing processes 
are much larger than time scale of microscopic particle interactions, as a result collapse can 
be divided into infinitesimal steps so that at each step the system is in equilibrium state 
which is given by Eq. (j4j). By using the chemical equilibrium Eq. (j4]) and conservation of 
total electric charge number, total baryon number and total lepton flavor number, number 
of each species of particles obtained as a function of collapsing star density p. We investigate 
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the collapsing system from white dwarf density p w d = 8.73 x 10~ 8 p nuc to the maximal density 
at the Schwarzschild radius of the star i.e., R s = 2GM, at which black hole is formed for 
stars with the mass M > 3.2 M . Maximal density p s = M/(^irR^) is plotted as a function 
of the mass M of collapsing stellar cores in Fig. H] to compare with the density thresholds 
p M = 2.7p nuc and p^ = 4.6p nuc of muon and pion productions. Based on this figure, there 
are three cases of collapsing stars. Case 1: stars with the mass 3.2 M < M < 3.8 M & can 
reach to the density thresholds p M and p n of both muon and pion productions. Case 2: stars 
with the mass 3.9 M < M < 4.9 M can only reach to the density threshold p^ of muon 
production. Case 3: stars with the mass M > 5M cannot reach to the density thresholds 
p M and p n of both muon and pion productions. According to these density thresholds there 
are four phases for neutrino production. The phase a begins from white dwarf density p„d to 
the transition density px = 2.8 x p nuc ( TTT]) . In this phase, via the inverse /3-decay reac- 
tion OH]) electrons and protons having been converted to the neutrons and electron-neutrinos. 
The Number of these neutrinos is given by Eq. ( )38l) . The phase b begins from the transition 
density px to the muon production density threshold p M for cases 1 and 2, and for case 3 
to the maximal density p s . In this phase, via /3-decay reaction (E]) neutrons are converted 
into protons, electrons and electron-antineutrinos. The number of these antineutrinos is 
shown in Fig. dSJ). The phase c begins from the muon production density threshold p M to 
the pion production density threshold p n for case 1 and for case 2 to the maximal density 
p s . In this phase, in addition to the /3-processes ((6]) and (0), through reactions (1T8]) and 
(II 9p . electron-antineutrinos and muon-antineutrinos are produced. The numbers of these 
antineutrinos are shown in Fig. ([6]). The phase d, which can only be reached in case 1, begins 
from the pion production density threshold p n to the maximal density p s . In this phase, 
Eqs. fl33|) and ( l34"j) show that by increasing the density p, electron number N e and muon 
number decrease. Based on lepton flavor number conservation, electron-neutrinos and 
muon-neutrinos are produced and their numbers are shown in Fig. [71 The ratio A/N p is an 
important quantity to show charge composition of compact star as a function of the density 
p. As muons and pions produced, this quantity decreases more rapidly as a function of the 
density p, as shown in Figs. [T]l3j Suppose that 10% of variation in gravitational energy fj49|) - 
(15T|) is converted to neutrino energy, we use the resulted number of neutrino productions (see 
Sec. IHIj) and Eq. ( 152]) to obtain the mean energy of neutrinos. Fig. |8] shows that the mean 
energy of produced electron-neutrinos in the phase a is about 10 MeV, which is in agreement 
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with neutrino energy observed in supernova explosions. Figs. l9lfTTl respectively show that 
the mean energy of electron-antineutrinos produced in the phase b, electron-antineutrinos 
and muon-antineutrinos produced in the phase c and electron-neutrinos and muon-neutrinos 
produced in the phase d, are about GeV. The mean energy of total neutrinos and antineu- 
trinos produced in the phases a, b, c and d as shown in Fig. [12] is about 90 MeV. The ratio of 
the total muon neutrinos to the total electron neutrinos M e i.e., = N^/N^ have been 
computed. For the flux of 10 MeV electron-neutrinos, which are produced in the phase a, 
and collapsing stars with the mass M > 5 M & this quantity is zero (0 = 0) and by con- 
sidering the effects of neutrino oscillation becomes 0( D - ) = 0.39 at detectors on the Earth. 
And for the flux of GeV neutrinos, which are produced in phases b, c and d, the quantity 
(0 ( - D - ) ) is computed at the source (at detectors on the Earth), the Fig. [TBI shows that < 1 
(0(°) < 1) at the source (at detectors on the Earth). This result, is completely different from 
high-energy neutrino production in compact stars without considering star as a completely 
degenerate Fermi gas which predict the ratio = 2, see processes (Op), at the source and 
0( D ) = 1 after oscillation at the Earth detectors. The number of antineutrino events in an 
ordinary detector such as Kamiokande is computed for a collapsing star at distance of 1 Mpc 
from the Earth. The number of GeV antineutrino events N cvcnt > 1, depending on the mass 
of collapsing stars, as shown in Fig. [14] with and without neutrino oscillation effects. Total 
energy of neutrino flux e > 10 53 erg, as shown in Fig. [15] 

We have to point out that in this approach it is not possible to determine produc- 
tion of neutrinos and antineutrinos via the neutral current interaction, for example NN- 
bremsstrahlung reaction NN — > NN + vv 33] . In addition, an other important assumption 
in this approach is that as soon as produced, all neutrinos escape away from collapsing stars. 
By using this assumption, we are able to approximately solve algebraic equations of equilib- 
rium conditions and conservations of particle numbers. In fact, some neutrinos are trapped 
and participate interactions inside collapsing star before they escape away from collapsing 
stars, for instance, reactions (fTBl) and ( 1301) . This assumption is on the basis of the chem- 
ical potential of trapped neutrinos being smaller than the chemical potential of electrons 
and muons. If the effects of neutrinos trapped are considered, the number and energy of 
neutrino emission should be smaller than estimates given in this article. Our results should 
be considered as approximate estimations and further numerical calculations are necessary. 
Nevertheless, the approximate results presented in this article give some physical insight 
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into the neutrino emission in gravitational collapses and show number- and energy-fluxes of 
neutrinos that could be relevant to observations. 
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